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ANALOG COMPUTER STUDY OF DESIGN PARAMETER EFFECTS ON THE 


STABILITY OF A DIRECT-ACTING GAS PRESSURE REGULATOR 

by Miles O. Dustin 
Lewis Research Center 

SUMMARY 

An analog computer study of a direct-acting gas pressure regulator is presented. 

The effect of regulator design parameters on stability, transient response, and steady- 
state accuracy were examined. The regulator used in the study is a single-stage, 
spring-loaded regulator used in the gas management system of a solar Brayton cycle 
power generator. The program was carried out in three parts: 

1. An investigation of the effect of poppet stem friction, sensing orifice diameter, 

and output volume on the dynamic performance of a regulator with a linear 
spring and simple poppet 

2. An investigation of the effect of an orifice in series with the poppet on regulator 

performance 

3. An analog simulation of the Brayton cycle regulator which includes a nonlinear 

spring and orifice in series with the poppet and comparison of the simulation 
with performance curves taken from tests performed on the actual regulator 

Poppet flow forces were determined experimentally. Flow forces were measured 
for both the simple poppet and for the poppet with a series orifice. It was found that the 
simple poppet flow forces work in a direction to aid in opening the poppet as the poppet 
is opened. This tends to make the regulator unstable. By adding an orifice in series 
with the poppet this effect is reversed. The force acts in a direction to close the poppet 
as it is opened, which has a stabilizing influence on the regulator. The orifice degrades 
steady-state regulation, however. 

It was also found that with no friction on the poppet stem, a large spring rate was 
required to stabilize the regulator. The large spring rate causes poor steady-state reg- 
ulation. However, if a small amount of friction is applied to the stem, the spring rate 
could be greatly reduced, resulting in improved regulation. Other factors which improve 
stability without affecting steady-state regulation are a decreased sensing orifice diam- 
eter and a larger output volume. 



INTRODUCTION 


The gas pressure regulator considered in this study was based on a regulator used 
in the gas management system in a solar Brayton cycle power generator. Since stabil- 
ity problems were encountered during the first checkout of the regulator, a study was 
carried out to determine those design parameters that affect the regulator stability and 
to determine the effect of these parameters on transient performance and steady-state 
regulation characteristics. 

The simple, direct-acting, gas pressure regulator is used in almost every pneu- 
matic system in use today. The modern industrial gas pressure regulator, having 
evolved over many years of application experience, provides good accuracy and stability 
characteristics over the range of pressures and flows of most commercial applications. 
However, the application of a given design to new pressure and flow conditions is uncer- 
tain and can result in unsatisfactory regulation. 

The transient behavior of the simple hydraulic pressure regulator was analyzed by 
Gold and Otto in reference 1 using an electrical analogy technique. The steady-state 
characteristics of the simple gas regulator have been studied by Iberall in reference 2. 
Also, a dynamic analysis of the simple regulator has been done by Tsai and Cassidy in 
reference 3. Both linear and nonlinear problems were considered in the latter. 

The difficulties encountered in making a dynamic analysis of the gas regulator are 
largely due to the nonlinearities of the equations describing the regulator. The nonline- 
arities include orifice flow through the metering valve, friction on the poppet stem, and 
the valve flow area as a function of stem position. The mathematical solution of the non- 
linear equations describing the dynamic behavior is difficult and time consuming. It is 
usually necessary to linearize the equations about an operating point and to assume only 
small deviations of the variable from this operating point. It is therefore, not possible 
to study the transient behavior of the regulator to large changes in output flow. 

The analog computer can be used to great advantage in this type of investigation. 

The nonlinear equations describing the regulator can be retained. Another important 
advantage of the analog computer is the ability to rapidly change the design parameters 
and to assess the effect of these changes on the dynamic performance of the regulator. 
This report describes an analog computer simulation of a simple, direct-acting gas 
regulator. The simulation includes the nonlinearities encountered in the real regu- 
lator. 

The investigation was conducted in three parts: 

(1) The first part, referred to as Case I in this report, considered the Brayton cycle 
regulator with two simplifications. A linear spring and a simple poppet metering valve 
were used. Changes were made in the design parameters to determine the effect of the 
changes on stability, transient performance, and steady-state regulation. 
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(2) To obtain good stability characteristics in the actual Brayton cycle gas regulator, 
the manufacturer added an orifice in series with the metering poppet valve. The second 
part of the study, or Case II, was to determine the effect of this orifice on the transient 
response and steady-state regulation of the regulator. All other fixed parameters such 
as bellows area, internal volumes, and valve area were the same as in the Case I study. 
As in the Case I study, changes were made in the design parameters to determine their 
effect on the regulator performance. 

(3) The third part, Case m, was a simulation of the Brayton cycle gas regulator, 
with a nonlinear spring and the orifice in series with the poppet valve. The performance 
of the simulated regulator was then compared to performance tests made on the actual 
regulator by the regulator manufacturer. 

REGULATOR DESCRIPTION 
General Description 

The simple, direct-acting gas regulator used for this study is shown schematically 
in figure 1. A bellows senses the regulated pressure P r . (All symbols are defined in 
appendix A.) The bellows force, acting against a spring, adjusts the flow through the 
poppet valve to maintain the regulated pressure at the desired value. The desired pres- 
sure is set by the initial load on the spring. 

This regulator has a balanced poppet valve. The poppet has been made insensitive 

to the regulators supply pressure P by the addition of a small bellows attached to the 

s 

bottom of the poppet. The regulated pressure P r is supplied to this bellows by means 
of a bleed hole drilled through the poppet valve stem. This allows the regulated pres- 
sure to act on both ends of the poppet with equal areas, thus balancing the poppet. 

Figure 2 more clearly illustrates the operating principle of the actual regulator. In 
this figure the pressure sensing bellows has been replaced by a piston. Pressure acts 
on areas Aj, A 2 , and A g - A^. These areas, which are separated by close clearances 
in figure 1, are separated by orifices A^ 2 and A^ in figure 2. The balancing bellows 
in figure 1 has been replaced by a piston which has the same diameter as the valve seat. 
The series orifice used in the Case II and Case III studies is shown by dashed lines. 

If only the regulated pressure and the initial spring force acted on the poppet stem, 
ideal pressure regulation would result. However, several other forces, which are func- 
tions of valve position and supply pressure, act on the system. These forces affect the 
steady- state accuracy of the regulator. These forces are listed and described below: 

(1) An additional force is furnished by the spring, due to stroking the poppet. As 
the flow rate or supply pressure to the regulator changes, a new poppet position is 
required. Since the spring load must assume a new value, the regulated pressure must 
also change. This change in regulated pressure contributed to the steady-state offset of 
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the regulator. The magnitude of the force is determined by the poppet displacement and 
by the spring rate of the spring. 

(2) Although the poppet is balanced, in that the same pressure is applied to both the 
top and bottom of the poppet, when flow occurs through the valve there is an increase in 
force on the top face of the poppet. This effect can be seen by referring to the two poppet 
conditons shown in the following sketch: 


P r 



Closed 



In the closed position the pressure P acts on the same area A„ on the top and bottom 
of the poppet. However, as the poppet opens, the location of the minimum flow area 
moves to a smaller diameter corresponding to A , as shown. Also, velocity through 
the poppet changes the static-pressure profile and allows P to act on a larger area. 

The net effect is to increase the force in a direction to open the poppet. (Additional in- 
formation on poppet valve characteristics can be found in ref. 4.) The flow force is ap- 
proximately proportional to valve stroke, just as is the spring force. However, the 
change in flow force is in the opposite direction, or in effect, is a negative spring rate. 
The spring rate of the spring has a stabilizing influence on the regulator, while the effec- 
tive negative spring rate of the flow forces can contribute to instability in the regulator. 
However, the change in flow force does help to compensate for the steady- state offset 
due to the change in spring force when the valve assumes a new steady-state position. 

The flow force will cause the poppet stem to lose contact with the pressure sensing bel- 
lows and to open fully unless the stem is attached to the bellows (fig. 1). Another spring 
inside the balancing bellows can also be used to force the valve stem to remain in contact 
with the pressure sensing bellows. 

(3) Sliding frictional forces on the valve stem resulting from guides on the stem can 
cause steady-state pressure errors. The frictional forces will also influence the stabil- 
ity characteristics of the regulator. 
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System Equations 


The general equations developed in this section were used in all three parts of the 
study. Changes to the equations necessary to perform the various studies were made as 
required. Figure 3 defines the bellows areas used in the equations. 

Dynamic equation . - The basic equation that describes the metering valve and bel- 
lows assembly dynamics is obtained by equating all forces acting on the assembly to zero. 
Positive forces are assumed to act upwards, and positive motion is assumed to be in the 
direction of opening the valve or downwards. The equation is 

Mx + K g x + PjAj + PgAg + P 3 ( A 3 " A 4 ) - F s “ + F f( s S n x ) = 0 (1) 

where 
Mx 


P 1 A 1’ P 2 A 2’ 
P 3 ( A 3 - A 4) 

F s 

Kjx 

F f (sgn x) 


+ F^ 


Velocity 



force due to accelerating mass M (mass M is the mass of all moving 
parts. ) 

additional spring force due to stroking metering poppet through dis- 
tance x 

pressure forces acting on bellows areas 
initial force set on spring 

flow force on poppet valve (For this program flow force constant Kj was 
obtained experimentally for each supply pressure value. Appendix B 
describes the experimental procedure for obtaining the coefficients.) 

friction force on valve stem caused by valve stem support guides (It is 
assumed that the friction force is of constant value, that it acts upward 
when the stem motion is downward, and that it acts downward when the 
stem motion is upward, as shown in the following sketch.) 
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Continuity equation for main regulator flow. - If it is assumed that flow into volume 
V is negligible, the continuity equation for flow into volume V r may be written as 


dP r 

dt 


= m - in 
v o 


( 2 ) 


where 



storage term for volume V r (C r is the pneumatic capacitance of volume V r 
and is equal to V r /kRT for an isentropic process. The effect of poppet stroke 
x on the capacitances C , Cj, C 2 , and C 3 is small and is neglected through- 
out the report.) 


flow through metering poppet (It is assumed that flow through the valve is 
always choked. ) For air, 


m v=' 


0. 532 C f P g A v 


/T 


(3) 


m v 


0.0400 C f P A. 


]/T 


(for use with SI units) 


(3a) 


where C f is the valve flow coefficient and is assumed to be 0. 8 (ref. 3), and 
A v is a function of poppet displacement. 

m^ regulator output flow which, since the output is always choked, is 


A o 


0. 532 C f P r A 0 


(4) 


m o 


0.0400 C f P r A 0 


\Tv 


(for use with SI units) 


(4a) 


Continuity equation for flow into volume Vj. - The continuity equation for the flow 
into volume Vj can be expressed as 


dP^ 

dt 


lij - m 2 + PjAjX 


(5) 
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where 

dp i 

C. storage term for volume V, 

1 dt 1 

mj flow into Vj from volume V r through orifice A g (This flow is assumed to be 

incompressible, which can be said if the pressure drop between volume Vj 

and volume V is small.) 

r ' 

* 1 = A s C f |/2 go P r |(P r - Pl )| [sgn(P r - P,i] (6) 

“2 ' A f2 C f Ao p r!< P l - P 2>l [ s e”( p l - P 2>] U> 

p r AjX equivalent flow rate into volume because of the bellows motion 

It is assumed that the densities p.p pg, Pg, and p f are equal. 

Continuity equation for flow into volume V2. - In a similar manner, the continuity 
equation for flow into V2 may be written as follows: 

dPn 

C 2 = m 2 - m 3 + P r A 2 x (8) 


where 


”3 = A f3 C f/ 2 <SoPrl< P 2 - P 3>| [ S «”< P 2 ' p 3»] < 9 > 

Continuity equation for flow into volume V3. - In a similar manner, the continuity 
equation for flow into volume Vg may be written as follows: 

dPo . 

Co = m 3 +P r (Ao - A 4 )x (10) 

dt 

Metering valve flow area . - The change in the metering valve flow area A y with re- 
spect to valve stroke x is dependent on the valve shape. The shape used in this program 
was that of the Brayton cycle regulator valve and is shown in figure 4 . The flow area is a 
discontinuous function of stroke x, and may be expressed in general terms as follows: 
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When 0 < x< (r - r )/cos 9 sin 9 

■T* V 


A =7 r x sin 6 (2 r - x cos 9 sin 9 ) 

V 


(ID 


When (r p - r y )/cos 9 sin 9 < x< x max 


A v = 7r(r p + r V> 


fr x .v r v>] 

2 + (r - r v ) 2 ) 

L tan 9 


< 

J 


1/2 


(12) 


ANALOG SIMULATION OF REGULATOR WITH SIMPLE POPPET - CASE I 


Discussion 

The first case studied was that of the regulator with a linear spring and a simple 
poppet metering valve. In terms of the system dynamic equation (eq. (1)), this simply 
means that the spring rate K g is constant. The flow force coefficient Kj for the simple 
poppet was determined experimentally, as described in appendix B. Also, the sign of 
KfX is negative, as shown. 

The simulation was performed on a general-purpose electronic analog computer. The 
specific equations used to construct the simulation are developed in appendix C. 

The actual gas pressure regulator used in the Brayton cycle gas management system 
is required to regulate the system pressure under two output flow conditions. The trans- 
ient responses conducted on the analog simulation were, therefore, made by switching the 
output flow orifices between areas of 0. 10 and 0.004 square inch (0.65 and 0.026 cm^). 
This is roughly equivalent to the acceptance tests performed on the actual regulators. Al- 
though the regulator was designed for ultimate use with a mixture of helium and xenon, the 
preliminary tests were conducted with nitrogen. The simulation was done with 530 R 
(290 K) air, which has nearly the same fluid dynamic properties as nitrogen at the same 
temperature. 


Simulation Study Results 

The study was designed to determine the effect of various system parameters on the 
stability of the regulator. The parameters involved were the sensing orifice diameter 
D g , poppet stem friction F f , output volume V r , and supply pressure P g . Table I gives 
a summary of conditions under which the study was conducted. The table also lists fig- 
ures in which the results appear. 
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Effect of stem friction. - The effect of various amounts of friction on the regulator 
stability was determined in the following manner: With a fixed supply pressure P g , 
sensing orifice diameter D g , output volume V r , and stem friction F^, the system was 
made unstable by lowering the spring rate. The spring rate was then increased until the 
system was stable. It should be noted that the system could always be made stable with 
a sufficiently high spring rate. The system was defined as stable if oscillations were 
not sustained following a step change in output flow orifice area. 

Figure 5 shows how the value of spring rate required to stabilize the system varies 
as a function of the stem friction F^. From these curves it appears that small amounts 
of friction on the poppet stem tend to stabilize the regulator. This is particularly true 
when the sensing orifice is large. When the sensing orifice is small, the system is then 
insensitive to small amounts of friction but can be unstable if the friction becoms large. 

Effect of sensing orifice diameter . - Figure 6 is a cross-plot of the data of figure 5 
and shows how the amount of spring rate required to maintain a stable system varies as 
a function of sensing orifice diameter for various values of friction. These curves also 
show that, if the sensing orifice is large, the system can tolerate large amounts of stem 
friction. However, for orifice diameters larger than 0.070 inch (0. 18 cm), at least a 
small amount of friction is very desirable to maintain stability. 

Transient responses for various values of friction and sensing orifice diameter are 
shown in figures 7(a) to (e). These responses were all conducted with the same spring 
rate, supply pressure, and output volume. In figure 7(a) it can be seen that a system 
which is initially unstable can be stabilized by decreasing the size of the sensing orifice. 
However, if the orifice is made too small, large overshoots in the transient response can 
result. 

Effect of output volume. - As the output volume is made larger, the system becomes 
more stable and can tolerate a lower spring rate. This effect is shown in figure 8, which 
is a plot of the spring rate necessary to stabilize the system as a function of the output 
volume. The largest sensing orifice diameter and smallest stem friction values were 
used to minimize their damping effects. It is desirable to operate the regulator with a 
minimum system spring rate, since this results in minimum steady-state regulation off- 
set. Figure 9 shows the transient response of the system to step changes in output flow 
orifice area for various output volumes. 

Effect of supply pressure . - With a linear spring, the regulator’s ability to maintain 
a constant regulated pressure is greatly affected by supply pressure. This effect is 
shown in figure 10, which shows the transient response for step changes in output flow 
for several supply pressures. Since the gain of the poppet valve is a direct function of 
supply pressure, the regulator with a linear spring is more stable at low supply pres- 
sures. This means that the value of spring rate must be selected to give good stability at 
the highest supply pressure. This results in poor regulation or greater offset at low sup- 
ply pressures, as shown in the figure. The ideal system should have a nonlinear spring 
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rate in which the rate changes as a function of stem position. When the poppet is nearly 
closed, which is the case for high supply pressure, the rate would be high. When the 
supply pressure is low and the poppet valve is further open, the spring rate would be 
lower. As will be shown later, this was done in the actual gas regulator considered in 
Case m. 


ANALOG SIMULATION OF REGULATOR WITH ORIFICE IN SERIES 
WITH POPPET VALVE - CASE H 
Discussion 

During the initial checkout of the Brayton cycle gas pressure regulators, the manu- 
facturer found that the regulators could be made stable by adding an orifice in series with 

the poppet valve. The orifice was placed downstream of the poppet valve, as shown in the 

o 

inset of figure 1. The orifice area was 0.041 square inch (0.26 cm ). This part of the 
analog study was to investigate the effect of the series orifice on regulator stability, tran- 
sient response, and steady-state regulation. The effect of adding the orifice can be seen 
by referring to the inset of figure 1. The orifice produces a restriction in the normal 
poppet flow such that the pressure on top of the poppet face varies as a function of poppet 
stroke if the regulated pressure P r is constant. The pressure on top of the poppet is 
fed to the bottom of the poppet and acts in an upward direction on the poppet. However, 
the area under the poppet is greater than the area on top of the poppet by an amount equal 
to the poppet stem. Thus, a net force is created that increases in an upward direction as 
the poppet is opened. The effect is in the same direction as the incremental spring force 
on the poppet and is in a direction to improve the regulator stability. The value of the 
flow force coefficient was determined experimentally, as described in appendix B. In 
equation (1), the addition of the series orifice causes the sign on the K^x term to be 
positive. 


Simulation Study Results 

The study to determine the effect of the orifice in series with the poppet was con- 
ducted on the analog computer in a manner similar to the Case I tests wherein the simple 

poppet was used. That is, the output orifice areas were switched back and forth between 

o 

areas of 0. 10 and 0.004 square inch (0.65 and 0.026 cm ). (A summary of conditions 
under which the study was conducted is given in table I. ) 
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Effect of supply pre ssure and friction . - An important difference in the behavior of 
the regulator with the series orifice was that the regulator was nearly always stable, 
even with very low spring rates. Also, since the value of the flow force coefficient 
increases with the supply pressure, the regulator is now more stable at high supply 
pressures. This effect can be seen by comparing the transient response curves shown in 
figures 11(a) to (c), which show transient behavior of the simulated regulator for various 
supply pressures and amounts of stem friction. The effect of stem friction on the per- 
formance of the regulator is also shown. If the regulator is quite unstable with no fric- 
tion, as shown in the top trace of figure 11(a), additional friction increases the amplitude 
of the oscillations. However, if the regulator is only slightly unstable, as shown in fig- 
ure 11(b), a small amount of friction stabilizes the regulator. If additional friction is 
present, the regulator output can assume additional error within limits determined by 
the friction value. The band of this additional error will be equal to the frictional force 

divided by the effective area of the bellows. The effective area of the bellows used in 

2 

this study was 2.0 square inches (13 cm ). Therefore, the output pressure can be ±5.0 
2 

psi (3.4 N/cm ) from the zero friction output value if the stem friction is 10 pounds 
force (44 N). 

Effect of sensing orifice diameter. - The effect of the sensing orifice diameter on 
the transient response to step changes in output flow is shown in figures 12(a) to (c). The 
effect is about the same as shown before for the simple poppet simulation. A system 
that is initially unstable can be stabilized by decreasing the sensing orifice size. How- 
ever, if the orifice is made too small, large overshoots can result. 

Effect of output volume. - Figures 13(a) to (c) show the effect of increasing the out- 
put volume. The effect is to make the system more stable. The system time constant, 
however, is increased. 


ANALOG SIMULATION OF BRAYTON CYCLE REGULATOR - CASE HI 

Discussion 

The first two parts of this study, Case I and Case II, evaluated on the analog com- 
puter the effect of many of the regulator design parameters on the stability, transient 
response, and regulation capabilities of the regulator. Many of these effects were also 
determined experimentally by the Brayton cycle gas regulator manufacturer during de- 
velopment and early checkout of the regulators. In addition to stabilizing the regulator 
by adding an orifice in series with the metering poppet valve and determining the optimum 
sensing orifice size, the regulator manufacturer installed a nonlinear spring. The effect 
of the nonlinear spring was to improve the low- supply-pressure offset by reducing the 



spring rate at large poppet openings. This was done in the actual regulator by using a 
Bellville-type spring which has nonlinear characteristics. The nonlinear spring charac- 
teristics shown in figure 14 were added to the analog simulation of the regulator by the 
addition of a diode function generator. The regulator manufacturer measured the spring 
force as a function of valve stroke in an assembled regulator. The value of stem friction 
was arbitrarily chosen as 0.20 pound force (0. 88 N) for the simulation. 


Simulation Study Results 

Figures 15 (a) to (c) show comparisons of simulated regulator performance with the 
actual regulator performance. The traces representing the actual regulator performance 
were redrawn from a light-beam oscillograph record. They were, therefore, expanded 
both in amplitude scale and time scale to match the analog simulation scales. They 
should, therefore, not be considered as accurate reproductions. Those features that are 
accurate are the magnitudes of the steady-state pressure levels for both high and low 
flows, the magnitudes of the pressure peaks, and the time for the pressure to reach the 
steady- state value following a transient. Only these factors will be considered in the 
discussion that follows. A summary of conditions under which the study was conducted 
is given in table I. 

Figure 15(a) compares the simulated regulator responses to step changes in output 
flow with similar responses of the actual regulator when the output volume was 100 
cubic inches (1600 cm ) for both. The simulated regulator curves follow the same gen- 
eral shape as the actual regulator. However, the actual regulator exhibited a higher os- 
cillation frequency. The overshoots when switching to low flow for the simulated regu- 
lator are higher than for the actual regulator. High overshoots were generally charac- 
teristic of small sensing orifice holes. This effect is demonstrated in figure 12(a). 

This could mean that the actual flow coefficient was greater than the 0. 7 value used in 
the simulation. The offset was greater for the simulated regulator than for the actual 
regulator. This means that the effective spring rate was higher in the simulation. The 
effective spring rate is a combination of the spring rate of the spring and the effect of 
the flow force on the poppet. 

O 

In figure 15(b) the output volume was increased from 100 cubic inches (1600 cm ) to 
400 cubic inches (6500 cm^) for both regulators. The overshoots are still larger in the 
simulated regulator, but it is not as noticeable with the larger output volume. Both reg- 
ulators are more stable, as expected. 

Figure 15(c) shows the effect of increasing the sensing orifice size from 0.025 inch 
(0.064 cm) to 0.036 inch (0.091 cm), while keeping the output volume at 400 cubic inches 
(6500 cm ). The simulated regulator has become more stable, while the actual regulator 
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has become less stable. As can be seen by referring to figure 6, the effect of changing 
the sensing orifice size depends upon the amount of stem friction. The value of stem 
friction for the simulation was chosen arbitrarily as 0. 2 pound force (0.88 N), which may 
not agree with the actual regulator. This means that changing the sensing orifice size 
may not affect the simulation in the same manner as it affects the actual regulator. 


SUMMARY OF RESULTS 

A single-stage, gas pressure regulator consisting of a spring-loaded, bellows- 
actuated poppet valve was analyzed. System equations were written and simulated on the 
analog computer. The computer simulation was used to determine the effects of the sig- 
nificant design parameters such as spring stiffness, stem friction, sensing orifice size, 
and output volume on the transient behavior of the regulator. Also represented in the 
simulation were certain design refinements such as a nonlinear spring, a poppet pressure 
balancing bellows and a downstream orifice between the poppet metering area and the 
pressure sensing orifice. The effects of the significant design parameters and refine- 
ments on regulator stability, response, and steady-state accuracy were demonstrated. 
Results from the analog computer simulation were compared with experimental results 
for a regulator designed for pressurizing a solar Brayton cycle space power generator. 
The mathematical model and computer circuit presented can be readily adapted to other 
similar regulator designs. Also, the trends pointed out could be used as guidelines by 
designers for making judicious improvements in the performance of regulator designs 
without recourse to detailed analysis. 

The following is a brief summary of the results of the study: 

1. An increase in system spring rate tends to improve stability; however, the regu- 
lating ability of the regulator is decreased. 

2. Decreasing the size of the sensing orifice tends to stabilize the regulator until an 
optimum orifice size is reached. Further reduction can cause instability, especially if 
the stem friction is great. It has no effect on the steady- state offset. 

3. Increasing the regulator's supply pressure can cause instability with a fixed 
spring rate. A nonlinear spring can provide good stability at all supply pressures with 
minimum steady-state offset. 

4. An orifice in series with the metering poppet valve can be used to improve the 
regulator's stability. The effect is similar to the nonlinear spring in that the orifice 
creates a stabilizing force which is greatest at high supply pressures. It also increases 
the steady- state offset. 

5. A small amount of friction can be of great value in stabilizing the regulator, es- 
pecially if the sensing orifice is too large to have a stabilizing effect. Too much friction, 
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however, can cause instability if the sensing orifice diameter is small. It can also cause 
the regulator to have a dead band whose width is a function of the friction value. 

6. Increasing the output volume has a definite stabilizing effect. The only detrimen- 
tal effect of increasing the volume is that the time constant of the system is increased. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 30, 1970, 

720-03. 
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APPENDIX A 


SYMBOLS 

2 2 

area, in. ; cm 

n o 

flow area between volume Vj and volume V 2 (see figs, land 2), in. ; cm 

flow area between volume V 2 and volume V 3 (see figs. 1 and 2), in. 2 ; cm 2 

2 2 

area of poppet measured at minimum flow area, in. ; cm 

2 2 

area of output orifice (see fig. 1), in. ; cm 

2 2 

area of poppet seat (see fig. 1), in. ; cm 

2 2 

area of sensing orifice (see fig. 1), in. ; cm 

2 2 

stem area used in poppet valve tests, in. ; cm 

2 2 

metering area of poppet (see fig. 1), in. ; cm 
capacitance, (in. 2 )(lbm)/lbf ; (cm 2 )(kg)/N 
flow coefficient 

diameter of sensing orifice (see fig. 1), in. ; cm 
friction force, lbf; N 
measured poppet force, lbf; N 
spring force, lbf; N 

gravitational constant, 386 (lbm)(in.)/(lbf)(sec 2 ); 100 (kg) (cm) /(N) (sec 2 ) 

flow force coefficient, lbf/in. ; N/cm 

spring rate, lbf/in.; N/cm 

ratio of specific heats 

mass, lbf-sec /in.; N-sec /cm 

mass flow rate (see fig. 1), lbm/sec; kg/sec 

mass flow rate from output orifice under high-flow -rate conditions (see fig. 4), 
lbm/sec; kg/sec 

mass flow rate from output orifice under low-flow-rate conditions (see fig. 4), 
lbm/sec; kg/sec 

mass flow rate through metering poppet (see fig. 1), lbm/sec; kg/sec 
mass flow rate through output orifice (see fig. 1), lbm/sec; kg/sec 



2 

P pressure, psia; N/cm 

R gas constant, in. -lbf/(lbm)(°R) ; cm-N/(kg)(K) 

Tp radius of poppet valve seat (see fig. 4), in. ; cm 

r y inside radius of bevel on poppet (see fig. 4), in. ; cm 

sgn[ ] sign of [ ] 

T temperature of gas, °R; K 

t time, sec 

3 3 

V volume, in. ; cm 

x stroke, in . ; cm 

x velocity, in. /sec; cm/sec 

x acceleration, in./sec^; cm/sec^ 

9 poppet bevel angle (see fig. 4), deg 

3 3 

p density, lbm/in. ; kg/cm 

Subscripts: 

b pressure balancing bellows 

f flow 

max maximum 

r conditions in output volume (see fig. 1) 

s supply 

1, 2, 3, 4 conditions under various bellows areas (see figs. 1 and 2) 
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APPENQfX B 


DETERMINATION OF FLOW FORCE COEFFICIENTS 

The flow force appearing in equation (1) as KjX, results from increased pressure 
loading on the top of the poppet as it is opened. In the closed position, equal pressures 
act on the top and bottom of the poppet. However, as the valve moves open, the mini- 
mum flow area moves to a smaller diameter. This allows the supply pressure to act on 
a larger area and the regulated pressure to act on a smaller area, resulting in a net in- 
crease in downward force on the poppet. Reference 4 gives emperical flow force equa- 
tions derived from experimental data of reference 3 for several poppet designs. In this 
study, however, it was decided to obtain the flow forces experimentally because of the 
unusual poppet design involved. 

A special test fixture, shown in figure 16, was built for obtaining the flow forces for 
both the simple poppet and for the poppet in series with an orifice. The orifice sleeve 
was soldered into the valve seat after the tests on the simple poppet configuration were 
completed. 

The poppet opening was varied by means of an adjustment screw on the force cell. 
The regulated pressure was controlled by means of a pneumatic operated valve on the 
output line. The regulated pressure was measured with a pressure transducer and the 
pressure was regulated by a proportional-plus- integral electronic controller. The out- 
put of the controller was converted to a pneumatic signal by means of a torque-motor- 
actuated, three-way valve. 

Tests were conducted in the following manner: The main valve from the nitrogen 
bottles was initially closed and the six valves on the nitrogen bottles were opened. The 
poppet opening was adjusted to the desired value. The main valve was then opened sud- 
denly, discharging the nitrogen bottles through the poppet valve fixture. Variables re- 
corded on an eight-channel strip-chart recorder include (1) supply pressure P g ; 

(2) pressure which normally would be in the pressure balancing bellows P^, (3) regu- 
lated pressure P , and (4) force on the load cell F m . 

To convert the force measured during these tests to an equivalent force for use in 
the simulation, it was necessary to add to the measured force F m a force equal to the 
stem area A g ^. times the regulated pressure P r . It is also necessary to subtract a 
force equal to the effective area of the balancing bellows times the measured bellows 
pressure P^. Figure 17(a) shows the resulting flow force plotted as a function of poppet 
stroke for the simple poppet. Figure 17(b) shows the flow force as a function of poppet 
stroke for the poppet in series with the orifice. Values of flow force coefficient 
were taken as the slopes of these curves at the operating poppet opening value for the 
maximum regulator flow condition. The curves shown as dashed lines were extrapolated 
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from data taken at lower pressures . Even though the nitrogen bottles were initially 

9 

charged at 2000 psig (1380 N/cm ), the supply pressure at the fixture under flowing con- 

2 

ditions rarely exceeded 1200 psig (830 N/cm ) because of pressure losses through the 
valves on the bottles. 
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APPENDIX C 


DEVELOPMENT OF SPECIFIC SYSTEM EQUATIONS FOR ANALOG SIMULATION 

The general equations used to describe the gas pressure regulator are given in the 
main text in the section System Equations. The specific equations used to construct the 
analog computer wiring diagram, figure 18, will be developed in this appendix. The 
equations appearing in brackets are for use with the SI units as shown in the symbol 
list (appendix A). Potentiometer settings are listed in table II. 

Dynamic Equation 

The general equation describing the dynamics of the metering valve and bellows 
assembly is given by equation (1) : 

Mx + K g x + P-jAj + P 2 A 2 + P 3^ A 3 " A 4^ " F s “ Kf x + F f (sgn x) = 0 (1) 

The coefficients for this equation were obtained directly from the manufacturer’s detailed 
drawings for the regulator and are as follows: 

M = 1. 86x10” 3 lbf-sec 2 /in. (3.26xl0“ 3 N-sec 2 /cm) 

Aj = 0.770 in. 2 (4.96 cm 2 ) 

A 2 = 0.244 in. 2 (1. 57 cm 2 ) 

A 3 - A 4 = 0. 969 in. 2 (6.25 cm 2 ) 

F = 360 lbf (1600 N) 
s 

Fj = Variable 

K„ = Variable 
s 

Kj- is a function of x and P g and must be adjusted in accordance with experimental 

data for each P value. Using these values , equation (1) can be written 
s 
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-x — (0. 538xl0 3 K s )x + (0.414 x10 3 )P 1 + (0. 131x10 3 )P 2 + (0. 520x 10 3 )P 3 - 194xl0 3 


- (0. 538x10 3 )KjX + (0. 538xl0 3 )F f (sgn x) (Cl) 


x= (0.307xl0 3 K g )x + (1.52xl0 3 )P 1 + (0. 481xl0 3 )P 2 + (1. 92xl0 3 )P 3 - 490xl0 3 


- (0. 307xl0 3 )K f + (0. 307xl0 3 )F f (sgn x) 


(Cla) 


Continuity Equation for Main Regulator Flow 


The general equation for flow into volume V r is 


dP 


(2) 


Since C r = V r /kRT, V r = 122 cubic inches (2000 cm 3 ), and T = 530° R (290 K), C r = 

2. 56x10"^ (in. 2 )(lbm)/lbf [l. 68x10"^ (cm 2 )(kg)/N] . Also, for choked flow through the 
poppet and a flow coefficient of 0. 8, 


m v = 1.85xlO- 2 A y P s 


<C2) 




m v = l. 89X10- ^A v P s 


(C2a) 


For flow through the maximum output flow orifice area of 0.10 square inch (0.6a cm; 
with a flow coefficient of 0.6, equation (2) becomes 


dP 

dt 


- = 0.723xl0 2 A v P g - 5.42P r 


(C3) 


dP r 

— = 11. 2A P - 5.42P 
dt v s r _ 


(C3a) 


For flow through the minimum orifice area of 0.004 square inch (0.026 cm ), equation (2) 
becomes 
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I 


= 0. 723xl0 2 A P - 0. 217P 

dt v s r 


(C4) 


— - = 11. 2A P - 0. 217P 
dt v s r 


(C4a) 


When a change in output volume was investigated, equations (C3) and (C4) were al- 
tered to account for the change in C r> This simply means that all coefficients in the 
equations are changed inversely proportional to the change in output volume. 


Continuity Equation for Flow Into Volume Vj 
The general equation for flow into volume Vj is 


dP. 


dt 


= mj - mg + PjAjX 


( 5 ) 


Since Cj = V^/kRT, Vj = 1.69 cubic inches (27.7 cm\ and T = 530° R (290 K), Cj = 

3. 55xl0“ 6 (in. 2 )(lbm)/lbf [2. 35xlO" 6 (cm 2 )(kg)/N]. 

The values of densities p r , pp p 2 , and pg were assumed constant and were calcu- 
lated at the set regulated pressure of 180 psia (124 N/cm 2 ). The area A-, was 0.77 

2 1 
square inch (4.96 cm ). Substituting these values into equation (5) yields 


dP 


dt 


1 = (2.82xl0 5 )m 1 - (2.82xl0 5 )m 2 +(1.15xl0 2 )x 


(C5) 


dP 


dt 


- = (4.26xl0 5 )m 1 - (4.26xl0 5 )m 2 + (0. 311xl0 2 )x 


(C5a) 


Flow mj into through the sensing orifice is assumed to be incompressible and 
may be obtained from equation (6) : 


“1 - A s C f Ac/>rl< p r - P l>| [ s g”< P r ‘ p l>] 


( 6 ) 
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Substituting in specific values of p r computed at 180 psia (124 N/cm ) 
value of Cj of 0. 7, 

m 1 = 0.45A s yj(P r - P^j] [f gn(P r - Pj)] 
m x = 3.8xl0" 2 A s y|(P r - P^j [sgn(P r - P^] 
Flow mg from volume Vj into Vg is obtained from equation (7) 
m 2 = A f gC f V 2 g 0 p r| (P l - P 2 )| [sgnCP! - Pgj] 


and an assumed 

(C6) 

(C6a) 

(7) 


2 

Substituting a flow area A^g of 0. 145 square inch (0.935 cm ) taken from manufacturers' 
drawings, equation (7) becomes 


m 2 = 0.06531/1^ - Pg)| [sgn^ - P 2 )] 
m 2 = S.SexlO^y^^Pj - Pg)| [sgn(P x - Pg)] 


(C7) 

(C7a) 


Continuity Equation for Flow Into Volume Vg 
The continuity equation for flow into volume Vg is 

dp 2 • • 

Cg = mg-m 3 +p r AgX (8) 


3 

Volume Vg measures 0.0127 cubic inch (0.208 cm ), which results in a value for Cg 
of 2.67X10" 8 (in. 2 )(lbm)/lbf [l.75xl0“ 8 (cm 2 )(kg)/N]. Substituting this value of C« and 

O “ 

a value for Ag of 0.244 square inch (1. 57 cm ) into equation (8) results in 


dPc 


dt 

dt 


(3. 75xl0 7 )mg - (3.75xl0 7 )rii 3 + (0. 484X10 4 ) x 
(5.70xl0 7 )m 2 - (5.70xl0 7 )m 3 + (0. 131xl0 4 )x 


(C8) 

(C8a) 
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Flow mg from volume V 2 into volume Vg is obtained from equation (9) : 


“3 - A f 3 C f V^o^r l^ P 2 " P 3^l f'S n ( p 2 P 3>] 


( 9 ) 


With a measured flow area A^g of 0.047 square inch (0. 303 cm j equation (9) becomes 


m 3 = 0. 0212 (P 2 - P 3 ) | [sgn(P 2 - Pg)] 


(C9) 


[™3 = 1.15xlO- 2 ^|(P 2 - P 3 )| [sgn(P 2 - P 3 )] 


(C9a) 


Continuity Equation for Flow Into Volume Vg 
The continuity equation for flow into volume Vg is 


dp 3 

C 3-T- = m 3 + M A 3- A 4) x 


( 10 ) 


O 

Volume Vg measures 0.112 cubic inch (1.84 cm ), which results in a value of C 2 of 

0. 235x10'® (in. 2 )(lbm)/lbf Jo. 155x10'® (cm 2 )(kg)/Nj. Substituting this value for Cg and 
the correct values for p r and Ag - A^ into equation ( 10 ) gives 


dPo c . o . 

“ = (4.26xlO b )m„ +(2.18xlO ,J )x (CIO) 

dt 6 


— = (6. 47xl0®)m„ +(5.93x10 2 x 
dt s 


(ClOa) 
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TABLE I. - SUMMARY OF CONDITIONS UNDER WHICH ANALOG COMPUTER STUDY WAS CONDUCTED FOR CASES I, H, AND III 


Case 


Objective 


Supply pressure, 
P s 


Spring rate, 
K 


Friction force. Output volume, Sensing orifice diameter. Figure 
F, V„ 

i r s 


psia N cm 


lbf in. N cm 


lbf 


Effect of friction force on stability | 2000 1380 

Effect of sensing orifice diameter on 
stability 

Effect of friction force and sensing 
orifice diameter on response 


Effect of output volume on stability 
Effect of output volume on response 
El feel of supply pressure on response Varied | Varied 


Varied 

Varied 

Varied 

Varied 

Varied 

Varied 

Varied 

Varied 

1235 

2160 

0 

0 


122 


2000 


.20 

.40 

2.0 

10 


.88 

1.78 

8.8 

44 


Varied 


Varied 


II 

EUoct ol friction force on response 

300 

207 

100 

175 

Varied 

Varied 

122 

2000 

0.125 

0.318 

11(a) 



1000 

690 





Varied 

Varied 





.125 

.318 

11(b) 



2000 

1380 





Varied 

Varied 





.125 

.318 

11(c) 


Effect of sensing orifice diameter 

300 

207 





0. 

20 

0.£ 

8 





Varied 

Varied 

12(a) 


on response 

1000 

690 













Varied 

Varied 

12(b) 



2000 

1380 









' 




Varied 

Varied 

12(c) 


El led of output volume on response 

300 

207 









Varied 

Varied 

0.125 

0.318 

13(a) 



1000 

690 









Varied 

Varied 

.125 

.318 

13(b) 



2000 

1380 









Varied 

Varied 

.125 

.318 

13(c) 

III 

El led of supph pressure on response 

Varied 

Varied 

Nonlinear 

Nonlinear 

0. 

20 

0.88 

100 

1600 

0.025 

0.064 

15(a) 



Varied 

Varied 

Nonlinear 

Nonlinear 

.20 

.88 

400 

6500 

.025 

.064 

15(b) 



Varied 

Varied 

Nonlinear 

Nonlinear 


20 

.88 

400 

6500 

.036 

.091 

15(c) 


5 

6 

7(a) 

7(b) 

7(c) 

7(d) 

7(e) 


T 

Varied 

T 

Varied 

.20 

.88 

T 

Varied 

T 

Varied 

T 

0.125 

T 

0.318 

8 

1040 

1820 

.20 

.88 

Varied 

Varied 

.125 

.318 

9 

1560 

2730 

.20 

.88 

122 

2000 

.125 

.318 

10 



TABLE II. - POTENTIOMETER SETTINGS FOR GAS 


REGULATOR SIMULATION 


Potentiometer 

Setting 

Potentiometer 

Setting 

1 

0.100 

18 

0.282 

2 

. 131 

19 

.115 

3 

.414 

20 

.282 

4 

.520 

21 

.100 

5 

a . 538 x10' 3 K 

22 

.180 

5 

b .538 

23 

F f 




6 

.194 


100 + F f 

7 

8 

.538 

, 538xlO" 3 K f 
.010 

24 

F f 



9 


100 + F f 

10 

. 570 

25 

.197 

11 

.361xl0" 3 P g 

26 

.375 

12 

5. SSXlO- 4 !*!,/^ 

27 

.048 

13 

5. 55xl0~ 4 m^/ C r 
. 100 

28 

.375 

14 

29 

.018 

15 

. 180 

30 

.218 

16 

. 100 

31 

.426 

17 

33.9 

32 

. 180 

s 

33^ 

.202 


a Case I and Case II only. 
b Case III only. 
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gure 3. - Schematic of bellows assembly showing areas used in developing 
dynamic equation. 





0 .01 .02 .03 .04 .05 

Valve stroke, x, in. 



Valve stroke, x, cm 

Figure 4. - Metering valve flow area as function of valve stroke. 



Spring rate required to stabilize system, K s , HI cm 



0 2.4 6 8 10 

Poppet stem friction, F f , Ibf 


0 10 20 30 40 

Poppet stem friction, F f , N 

Figure 5. - Effect of poppet stem friction on regulator stability. Case I, 
simple poppet; supply pressure, 2000 psia (1380 N/cm 2 ); output 
volume, 122 cubic inches (2000 cm 3 !. 



,02 . 04 . 06 . 08 .10 .12 


Sensing orifice diameter, D $l in. 


.05 .10 .15 .20 .25 .30 

Sensing orifice diameter, D $ , cm 

Figure 6. - Effect of sensing orifice size on regulator stability. Case I, simple 
poppet; supply pressure, 2000 psia (1380 N/cm 2 ); output volume, 122 cubic 
inches (2000 cm 3 ). 




(b) Friction force, 0. 20 pound force (0. 88 N). 


Figure 7. - Simulated regulator response to step changes in output 
flow orifice area for various sensing orifice diameters. Case I, 
simple poppet; supply pressure, 2000 psia (1380 N/crrr); output 
volume, 122 cubic inches (2000 cm 3 b spring rate, 1235 pound 
force per inch (2160 N/cm). 
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Regulated pressure, P r , N/cm' 


130 p 
120 — 
no — 

130 1 — 
120 — 
cm no — 

130 p 
120 — 
no — 

130 j — 
120 — 
110 — 

130 — 
120 — 
110 — 




130| 

120 

no 


130 

120 

110 


130 
13 

I 120! 

S’ 

“ no 1 — 


130 

120 | 

no 



Figure 7. - Concluded. 
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Spring rate required to 
stabilize system, K s> N/cm 



2 3 4 5 6 7 SxlO 3 

Output volume, V r> cnr 

Figures. - Effect of output volume on regulator stability Case I, simple poppet; 
supply pressure, 2000 psia 1 1380 cnr^i; sensing orifice diameter, 0. 125 inch 
(0. 318 cm); stem friction, 0. 2 pound force i0. 8S N). 


CM 

E 


130 

120 

no 


130 

120 

110 


130 

120 

110 


130 

120 

110 



Figure 9. - Simulated regulator response to step changes in output flow 
orifice area for various values of output volume. Case I, simple 
poppet; supply pressure, 2000 psia {1380 N/cm); stem friction, 

0. 2 pound force (0. 88 N); sensing orifice diameter, 0. 125 inch 
(0.318 cm). 
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Figure 10. - Simulated regulator response to step changes in output flow 
orifice area for various values of supply pressure. Case I, simple 
poppet; spring rate, 1560 pounds force per inch (2730 N/cm 2 ); stem 
friction, 0. 2 pound force (0. 88 N); sensing orifice diameter, 0. 125 
inch (0. 318 cml; output volume, 122 cubic inches (2000 cm 3 ). 
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Regulated pressure, P r , N/cnV 



(a) Supply pressure, 300 psia (207 N/cm 2 ). (b) Supply pressure, 1000 psia (690 N/cm 2 ). 



(c) Supply pressure, 2000 psia (1380 N/cm 2 ), 


Figure 11. - Simulated regulator response to step changes in output flow orifice area for various values of stem friction. Case II, 
orifice in series with poppet; spring rate, 100 pounds force per inch (175 N/cm); sensing orifice diameter, 0. 125 inch (0.318 cm); 
output volume, 122 cubic inch (2000 cm^). 
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Regulated pressure, P r , N/cm' 
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Figure 12. - Simulated regulator response to step changes in output flow orifice area for various sensing orifice diameters. Case II - 
orifice in series with poppet; spring rate, 100 pounds force per inch (175 N/cm); stem friction, 0.20 pound force (a 88 N); output 
volume, 122 cubic inches (2000 cm3). 
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I 



Regulated pressure, P r , N/cm' 


Output volume, 



(a) Supply pressure, 300 psia (207 N/cm 2 ). (b) Supply pressure, 1000 psia (690 N/cm 2 ). 



Figure 13. - Simulated regulator response to step changes in output flow orifice area for various values of output volume. 
Case II, orifice in series with poppet; spring rate, 100 pounds force per inch (175 N/cm); stem friction, 0. 20 pound 
force (0. 88 N); sensing orifice diameter, 0. 125 inch (0. 318 cm). 
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Figure 16. - Test setup for measuring flow forces on poppet -type valve. 
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